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Epitaxial strain energies of epitaxial films and bulk superlattices are studied via first-principles total 
energy calculations using the local-density approximation. Anharmonic effects due to large lattice 
mismatch, beyond the reach of the harmonic elasticity theory, are found to be very important in 
Cu/Au (lattice mismatch 12%), Cu/Ag (12%) and Ni/Au (15%). We find that (001) is the elastically 
soft direction for biaxial expansion of Cu and Ni, but it is (201) for large biaxial compression of Cu, 
Ag, and Au. The stability of superlattices is discussed in terms of the coherency strain and interfacial 
energies. We find that in phase-separating systems such as Cu-Ag the superlattice formation energies 
decrease with superlattice period, and the interfacial energy is positive. Superlattices are formed 
easiest on (001) and hardest on (111) substrates. For ordering systems, such as Cu-Au and Ag-Au, 
the formation energy of superlattices increases with period, and interfacial energies are negative. 
These superlattices are formed easiest on (001) or (110) and hardest on (111) substrates. For Ni- 
Au we find a hybrid behavior: superlattices along (111) and (001) behave like in phase-separating 
systems, while for (110) they behave like in ordering systems. Finally, recent experimental results 
on epitaxial stabilization of disordered Ni-Au and Cu-Ag alloys, immiscible in the bulk form, are 
explained in terms of destabilization of the phase separated state due to lattice mismatch between 
the substrate and constituents. 



PACS numbers: 62.20.Dc, 68.60.-p, 81.10.Aj 

I. INTRODUCTION 

Recently, there has been much interest ^"^"^ in growth 
of epitaxial metal films and superlattices due to their 
unusual physical properties. The quality and struc- 
ture of these systems is of paramount importance 
for applications. Epitaxial monolayer and multilayer 
(up to 10 layers) formation has been observed for 
many metal/semiconductor and metal/metal combina- 
tions. Most metal/metal superlattices have been grown 
for elements in different crystal structures (e.g., fcc/bcc) 
and with considerable size mismatch (e.g., 10% for 
Cu/Nb^^^^^). Furthermore, elemental metals and alloys 
have been found to form epitaxially in structures which 
are unstable in bulk form.-'^*"^^ Recently, the topic of 
surface alloy formation in bulk immiscible systems has 
attracted considerable attention. ^■^^^^ These systems are 
usually strained due to film/substrate lattice mismatch. 
One would like to understand and predict the stability of 
these types of strained materials. In order to do so, one 
requires knowledge of two types of energies. The stability 
of epitaxial Ai-^Bx alloy films and strained ApBq super- 
lattices depends on (i) the energies of coherently strained 
constituents A and B, and (ii) the formation energy of 
Ai-xBx or ApBq itself. Regarding (i), previous theo- 
retical studies^^'*'^"^^ have described these energies us- 
ing harmonic models, but we are interested here in large 
strains for which the harmonic theory could break down. 
Thus, we develop a generalization of previous methods to 
treat the anharmonic epitaxial strain energies of the con- 



stituents. Regarding (ii), these energies depend on the 
configuration degrees of freedom of the epitaxial film, so 
their calculation requires statistical methods. ^^'^^ In the 
present paper we investigate items (i) and (ii) above us- 
ing accurate first-principles LDA calculations. 

As for (i), the constituent strain energy, we find that 
the harmonic strain theory,^*'^^ predicting a single, uni- 
versal relation for elastically soft directions, breaks down 
for sufficiently large substrate/film lattice mismatch. We 
find that under biaxial expansion, noble metals are soft 
along (001), but that under compression the soft direc- 
tion changes to (201). It is shown that the softness of 
(001) is a consequence of low bcc/fcc energy differences in 
noble metals, while the softness of (201) under compres- 
sive strain can be explained by loose packing of atoms in 
the {201} planes. Furthermore, the elastic strain energy 
as a function of direction exhibits qualitative shifts in the 
hard and soft strain directions, which cannot be guessed 
from the harmonic elasticity theory. For instance, we 
find that (110) becomes the hardest direction under bi- 
axial expansion, and (201) becomes the softest direction 
under biaxial compression, while the harmonic theory al- 
ways predicts either (111) as the hardest and (001) as 
the softest direction, or vice versa. 

Regarding (ii), the formation energy, we find that the 
anomalous elastic softness of the constituents along (001) 
and (201) leads to low constituent strain energy in super- 
lattices along these directions, which makes them more 
stable than superlattices along other G. For instance, in 
the size-mismatched systems Cu-Au, Cu-Ag, and Ni-Au, 
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AnBn supcrlatticcs along (001) arc the most stable for all 
periods n. Interfacial energies are found to be negative 
in Ag-Au and Cu-An (reflecting their bulk miscibility) , 
and positive in the phase separating systems Cu-Ag and 
Ni-Au. However, attraction between (110) interfaces in 
Ni-Au is very strong and favors short-period (n oc 2) 
superlattices over long-period superlattices with few in- 
terfaces. 

In the case of cpitaxially grown disordered alloys, wc 
find that the biaxial constraint on the phase separated 
constituents may stabilize the alloy with respect to phase 
separation. The stabilization effect is always greater 
on substrates oriented along elastically hard directions 
(i.e., with high constituent strain energy) like (111) than 
along soft directions like (001). For instance, on lattice- 
matched substrates, epitaxial Nio.sAuo.s alloys are stable 
at all temperatures, and Cuo.sAgo.s alloys are stable for 
T > 150 K if grown on a (111) substrate, although both 
these systems phase separate in bulk form or if grown on 
a (001) substrate. These predictions agree very well with 
recent experimental observations. ^^'^^ 



A^7(ai, a2, G) = mm [EY{ai,a2,c)] - E'^\aA). (1) 

In what follows, we are interested in the case where both 
the substrate and the unstrained bulk element A have the 
fee crystal lattice. Then ai and a2 are proportional to the 
equilibrium unstrained lattice vectors of fee A, a)Oi{A): 



a A, 



(2) 



where and are fee lattice parameters of the sub- 
strate and A^ correspondingly. The epitaxial strain en- 
ergy becomes a function of the substrate lattice constant 
and direction G only: 

Ai?7[(a,/a^)ai, (a,/a^)a2, G] = b.Ej{a,, G). (3) 

LDA calculations of A£J^P'(as,G) are described in 
Sec. III. 



B. Formation enthalpies of alloys and superlattices 



II. BULK AND EPITAXIAL STABILITY 
CRITERIA 

The stability of either free-standing or coherently 
strained alloys and superlattices requires specification of 
(i) epitaxial strain energies of pure constituents due to 
film/substrate lattice mismatch, (ii) formation enthalpies 
of disordered alloys (with respect to either strained or 
unstrained bulk constituents) and superlattices. In this 
section, we define these quantities and discuss the phys- 
ical situations where they should be used. 



A. Epitaxial strain energies of elemental constituents 

We start by considering (i) above, which is a common 
element to alloys and superlattices. Consider a film of 
pure element A coherently strained on a substrate ori- 
ented along direction G with surface unit cell vectors 
ai and a2, orthogonal to G. Wc assume that the film, 
being much thinner than the substrate, maintains co- 
herency with the substrate and plastically deforms to ac- 
comodate the lattice mismatch at the interface. This 
assumption is valid for films thinner than the critical 
thickness for the nucleation of misfit dislocations. Fur- 
thermore, we consider films which are thick enough so 
that the chemical interaction energy at the film/substrate 
interface and film/vacuum surface is negligibly small in 
comparison with the elastic deformation energy of the 
film. Under these assumptions, the epitaxial strain en- 
ergy AE'^\ai,a2,G) of film A is the strain energy of 
element A deformed in the growth plane to the unit cell 
vectors {ai,a2} of the substrate, and relaxed with re- 
spect to the out-of-plane vector c: 



Like the formation enthalpy of any ordered bulk com- 
pound, the formation enthalpy AHgl^^{pq,G) of an 
ApBq unstrained (bulk) superlattice is defined as the en- 
ergy gain or loss with respect to unstrained bulk con- 
stituents: 



AHll'\pq,G)=E'°\Aj,Bg,G) 



ET{aA) + ^E'£\aB) 



P + Q 



p + q 



(4) 



where qa is the equilibrium lattice constant of the un- 
strained bulk element A and E*^*'(aA) is the total en- 
ergy of A. This enthalpy characterizes the propensity 
to form superlattices with respect to the phase sepa- 
rated bulk constituents. If AH^l^^{pq,G) < 0, the 
unstrained superlattices are energetically favored over 
the phase separation, while the phase separated state 
is favored if AH^l^^pq,G) > 0. To be stable, free- 
standing bulk superlattices must satisfy stability crite- 
ria with respect to at least: (i) phase separation into 
unstrained bulk constituents and (ii) formation of a con- 
figurationally disordered bulk alloy. The bulk mixing en- 
thalpy, AH^^{Ai_^B^), of the alloy is given by: 



AH^^^{Ai_,B,) = E'^^'iA^.^B,) 
-[{l-x)ET{aA)+xEr{aB)], 



(5) 



where x = q/ {p-\-q) is the composition and E^°^{Ai-xBx) 
is the total energy per atom of the configurationally ran- 
dom alloy. 

If AH^^{A^_M < AH^l'^^iApBg) < 0, then 
both the superlattice and disordered alloy are stable 
with respect to phase separation, but the superlattice 
is unstable with respect to disordering. However, if 
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AH^l'^iApBg) < AH^}^{Ai_M < 0, then superlat- 
tices arc stable with respect to both phase separation and 
disordering, and it may be possible to grow them. 

The bulk formation enthalpy of a superlattice, 
AHg^^^{pq,G), can be separated into two components. 
To identify them, it is useful to first consider the infi- 
nite period superlattice limit p,q oo, where A/B in- 
terfacial interactions contribute a negligible amount of 
order 0{l/p). In this case, the bulk formation enthalpy 
of AooBoo superlattice is given by 



mm 

ai,a2 



oo,G) = AE'c%i^,G) (6) 
(1 - x)AET{a,,aL2,G) + xAE''/{auSL2,G)\ , 



where AE^^ is the epitaxial deformation energy of A, 
given by Eq. (1). We define this energy as the "con- 
stituent strain" (CS) to emphasize that in this limit 
the superlattice formation enthalpy depends only on its 
strained constituents. This is also the energy required to 
keep A and B coherent. 

For finite-period superlattices, the formation energy is 
determined not only by the elastic strain energy, but also 
by interactions between unlike atoms at A/B interfaces. 
We define this interfacial energy I{pq, G) as: 



AifsT^(M, G) - AH^l'^ipq ^ oo, G) = 



I{pq,G). 
(7) 



It is the total energy per layer of a single interface 
between infinite slabs of A and B oriented along G. 
7(oo) < signals that the interface is energetically fa- 
vored, while /(oo) > indicates that an isolated interface 
is not preferred, and long-period superlattices with fewer 
interfaces are usually more stable than the short-period 
ones (however, this simple argument is not always true, 
see the following discussion). 

For equiatomic (A)„/(B)„ superlattices Eq. (7) be- 
comes: 

AHl^'^in, G) = + AE^clix = 0.5, G). (8) 

For small n interfaces will interact with each other. We 
describe this process by the interface interaction energy 
6I{n,Gy. 



dl{n, G) = I{n, G) - I{n ^ oo,G). 



(9) 



Negative 5I{n,G) may favor short-period superlattices 

over long-period superlattices even if the interfacial en- 
ergy I(n oo,G) is positive. For this to happen it is 
necessary that 



dI{n,G) < -|/(n^ oo,G)| 



(10) 



In Sec. IV C we show that this unusual phenomenon oc- 
curs in Ni-Au. 



If a disordered alloy is grown epitaxially on a lattice- 
matched fee substrate, its stability with respect to phase 
separation is given by the epitaxial mixing enthalpy: 



SH^^^iA^.M = AH^^^{A,.M 
-(1 - x)AEf{as,G) - xAEf\a,,G), 



(11) 



where AE'^^^{as,G) is the epitaxial strain energy of 
Eq. (3), accounting for the fact that the phase-separated 
consituents must also be lattice-matched with the sub- 
strate. Due to the presence of those terms, disor- 
dered alloys may form epitaxially [5H'^^^{Ai-xBx) < 0] 
even if the corresponding bulk alloys phase separate 
[AH^^{Ai^xBx) > 0]. This situation is especially likely 
to occur for clastically hard directions G with large val- 
ues of AE^^g{as,G), for instance (111) and (110) (see 
Sec III B). 

The objective of this work is to calcu- 
late AHl^^ApB,) [Eq. (4)], AH^^{A,^M [Eq. (5)] 
and 6H^:^{Ai-xBx) [Eq. (11)] from first principles for 
Ag-Au, Cu-Ag, Cu-Au, and Ni-Au. This requires: 

(a) Epitaxial strain energies of pure constituents, 
AEXia^G) [Eq. (3)], for Ag, Au, Cu and Ni. This 
is described in Sec. III. 

(b) Equilibrium constituent strain energy AiJ^g 
[Eq. (6)] for Ag-Au, Cu-Ag, Cu-Au, and Ni-Au. This 
is described in Sec. IV A. ^ 

(c) The interfacial energy I{pq, G) of Eq. (8) re- 
quires AH'^^^^{ApBq,G) for arbitrary pq and G. 

AH^^{Ai_M and SH^}^{Ai_M require the total 
energy of a configurationally disordered solid solution. 
All these quantities are obtained from the mixed-space 
cluster expansion as described in Sec. IV B. 



III. ELEMENTAL EPITAXIAL FILMS 

A. Anharmonic epitcixial strain in thin films of pure 
elements: Analytic forms 

The epitaxial strain energy [Eq. (3)] of a film of ele- 
ment A (with an equilibrium fee lattice constant a a) on 
a fee substrate with lattice constant Og, oriented along 
direction G, is conveniently obtained in a two-step pro- 
cess considered by Hornstra and Bartels.'*'^ First, the fee 
crystal of bulk A is uniformly stretched (or compressed) 
to the lattice constant of the substrate Og. The energy 
change relative to free A is given by the hydrostatic bulk 
deformation energy AEj^^^{aa). In the second step, out- 
of-plane imit cell vector c of the film relaxes to satisfy 
Eq. (1). The change Ac = c — (as/a^)c° (where is 
the fee lattice vector of unstrained A), has components 
parallel [Acy] and perpendicular [Ac_l] to the growth di- 
rection G. The parallel component Ac|| changes the vol- 
ume of the unit cell and thus has a large effect on the total 
energy. In contrast, the so-called shear strain Ac_l shifts 
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planes orthogonal to G and docs not change the volume of 
the unit cell. Consequently, it has a much smaller effect 
on the total energy. Furthermore, this strain vanishes by 
symmetry for directions (001), (111) and (110), and the 
shear strain energy must have zero angular derivatives at 
these points. Therefore, we neglect the shear strain Acj_ 
also for low-symmetry directions. Bottomley and Fons'^^ 
have shown that this approximation introduces rather 
small errors in the harmonic epitaxial strain energies. 

Neglecting the shear strain Ac_l, the strain energy of 
element A is then a function of the direction G and two 
scalar variables, as and e|| = |Ac|||/as — 1. The epitaxial 

strain energy Ai?^P'(as, G) of Eq. (3) is the minimum of 
the strain energy with respect to at a fixed substrate 
lattice constant a^: 

A^7(a«, G) = min Wx\as,e\\ , G)] - E'f{aA). (12) 

The epitaxial strain energy AE'^^Us, G) is related to the 
epitaxial softening function^^'^* q{ag,G) by the relation: 



q{as,G) = 



(13) 



where AE^^^{as) is the hydrostatic deformation energy 
of fee A to the substrate lattice constant a^. The function 
Eq. (13) quantifies energy lowering due to the relaxation 
of c{A) in the second step of the deformation process 
considered above. 

The harmonic elasticity theory without the shear 
strain glvcs^^'"'''''^^ qha.rm{G) which depends on the growth 
direction G but not on the substrate lattice constant as'- 



3harm(G) = 1 — 



B 



(14) 



Cii+A (G) 

where B = |(Cii+2Ci2) is the bulk modulus, A = C44 — 
^(Cii — C-12) is the elastic anisotropy parameter, and 
7harm(G') is a geometric function of the spherical angles 
formed by G: 



(15) 



Ki is the Kubic harmonic of angular momentum I. The 
equilibrium value of the C|| / a ratio of the film is given by 

cl'^{as,G) = as{l + e\\) = ua - [2 - 3gharm(G)](a,5 - qa)- 



(16) 



Expansion of Y(a^,G) 




For the principle high-symmetry directions we have 



3.2 3.4 3.6 3.8 4.0 4.2 

o 

Substrate lattice parameter (A) 



7harm([001]) = 0, 7harm([110]) = 1, 7harm([lll]) = 



3' 

(17) 



FIG. 1. Expansion coefficients bi{ae) of Eq. (18) for Ag, 
Au, Cu, and Ni. 
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A parametric plot of 7 is presented in Ref. 56, which 
shows that the minimum of 7(G) is along (001) and the 
maximum - along (111). Therefore, depending on the 
sign of the elastic anisotropy A, f/harmCG) is either low- 
est for the (001) direction, and then (/harmCllH]) is the 
highest, or vice versa. Other directions always have in- 
termediate values of qhiirm{G). 

If anharmonic effects are important, q becomes a func- 
tion of the substrate lattice parameter a^. As we will 
show in Sec. IIIB, for deformations 2{as — a a) /{as -^o,a) 
of approximately 4%, the "exact" LDA q{as,G) exhibits 
appreciable dependence on the substrate lattice param- 
eter as and certain qualitative features cannot be repro- 
duced by the harmonic functional form Eqs. (14)-(15). 
Furthermore, sufficiently large epitaxial strains may take 
the lattice from the face-centered cubic (fee) structure 
into other low-energy structures [e.g., body-centered cu- 
bic (bcc) and body-centered tetragonal (bet)], causing 
anomalous softening of q{as, G) for these directions. Sec- 
tion IIIB shows that this indeed happens for (001) epi- 
taxial strain when > a,A- Therefore, Eqs. (14) (1.5) 
must be generalized to account for nonlinear effects be- 
yond the reach of the harmonic theory. This is achieved 
by replacing in Eq. (14) 7harm(G) by 7(as,G), where 

'max 

7(a.,G) = 7harm(G) + bi{a,)Ki{G) (18) 

(=0 

includes higher Kubic harmonics. For cubic systems I = 
0, 4, 6, 8, . . . The general expression for q is 



q{as,G) = 1 - 



B 



Gn+A7(a„G) 



(19) 



We have chosen this particular form for 7 since it guar- 
antees that all expansion coefficients tend to zero in the 
harmonic limit: 



lim bi{as) 

la— >aA 



0. 



(20) 



In summary, to calculate AE'^P^{as,G) of Eq. (3) we will 

use Eq. (12) to obtain it from LDA for a few substrate 
lattice parameters and along selected symmetry direc- 
tions G. We will also need to obtain the harmonic elastic 
constants Cu , C12 and C44. The calculated Ai?°P'(a,5, G) 
results are then fitted by the general Eqs. (13), (18) and 
(19). 



B. Anharmonic epitaxial strain of thin films of pure 
elements: LDA results 




3.2 3.4 3.6 3.8 4.0 4.2 
Substrate lattice parameter (A) 



We have calculated the epitaxial strain energy 
Ai?®P'(as,G) for Cu, Ni, Ag and Au along six prin- 
ciple directions (001), (111), (110), (113), (201) and 
(221). The local-density approximation^* (LDA), as 



FIG. 2. The calculated epitaxial softening functions 
q{as,G) for Cu, Ni, Ag and Au. Points represent the di- 
rectly calculated LDA values and lines show the fit using 
Eqs. (19)-(18). 
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FIG. 3. Epitaxial softening function q(as,G) for (a)-(c) Cu and (d)-(e) Au, at different values of the substrate lattice 
constant as- 



implemented by the linearized augmented plane wave 
(LAPW) mothod^^, was used to obtain the total ener- 
gies in Eqs. (12) and (13). q{as,G) was calculated from 
Eq. (13) and fitted with the functional form Eqs. (18) 
(19). The angular momentum cutoff in Eq. (18) was set 
to imax = 10, leaving five independent coefficients for 
each value of the substrate lattice parameter Og. This 
choice allows reproduction of the LDA values with a max- 
imum error of 0.04. The calculations have been done for 
biaxial compression (as < Ueq) of Au and Ag, for bi- 
axial expansion (ag > Acq) of Ni, and for both biaxial 
expansion and compression of Cu. The expansion coeffi- 
ciets bi{as), entering Eqs. (18), are shown in Fig. 1. At 
the equilibrium lattice constant a^q (vertical arrows in 
Fig. 1), where the harmonic formula Eq. (15) is exact, all 
bi are exactly zero. As deviates from acq, they change 
rapidly indicating the importance of anharmonic effects. 
In Cu and Ni for Ug > a^q, I = 6 term is as important as 
I = and I = 4 terms, contributions from I > 8 being an 
order of magnitude smaller. In Au for a^ < acq, 6o(as) 
and bs{as) are the dominating terms, while the behavior 
of Ag is mainly determined by 5o(fls) and bc,{as). Thus, in 
spite of broad similarities between the studied elements, 
they exhibit some interesting differences. 

Figure 2 shows the calculated LDA epitaxial softening 
functions q{as,G) of Eq. (13) for Cu, Ni, Ag and Au. 
There are important qualitative and quantitative differ- 



ences between gharm(G) given by the harmonic elasticity 
Eq. (14), and the anharmonic q{as,G) calculated from 
the LDA. First, all q{as, G) depend on the substrate lat- 
tice constant a^, while the harmonic q'harm(G) are inde- 
pendent of Gs- Figure 3 shows the directional depen- 
dence of q{as, G) for Cu and Au at a few values of a^: 
equilibrium lattice parameter of Cu (3.56 A), equilibrium 
lattice parameter of Au (4.04 A), and halfway between 
them (3.83 A). By construction, g at as = acq is given 
by the harmonic form Eqs. (14)-(15), shown for fee Au 
in Fig. 3(c) and fee Cu in Fig. 3(d). Epitaxial deforma- 
tion of Au with as < aeq makes the lobes along (111) 
much more pronounced than in the harmonic case. Fur- 
thermore, q for Au develops additional lobes along (001), 
which in the harmonic approximation is the softest di- 
rection. In contrast, q of Cu under biaxial expansion 
exhibits pronounced deepening of the (001) minima, but 
develops maxima along (110). 

Second, in the harmonic elasticity theory of Eq. (14) if 
(001) is the softest direction (smallest gharm), then (111) 
must be the hardest direction, and vice versa. Figure 2 
shows that this order does not hold for large deforma- 
tions: the hardest direction in Ni and Cu for a^q 
is (110), while the hardest directions in Ag and Au for 
as < acq are (111) and (001), both (110) and (201) being 
much softer than the former. 



6 



Epitaxial (100) strain energy of Cu 




o 

-41^ 
CS 

> 

a. 



e 



0.8 0.9 1.0 1.1 
Tetragonal shear y 




3.4 3.6 3.8 4.0 4.2 4.4 
Substrate lattice constant (A) 



FIG. 4. Contour plot of the two-dimensional energy surface £(7, V) for Cu. The continuous line shows the epita^xial path 
determined by Eq. (23), while the dashed line is the relation V = V{'y) obtained by minimizing E{-y,V) with respect to the 
volume y at a constant 7. The right panel shows the epitaxial strain energy as a function of the substrate lattice constant in 
comparison with the (much larger) bulk deformation energy AE^^^^{as). 



Third, Fig. 2 shows that q{as, G) of different directions 



cross for substratc/filni lattice mismatch 2|as 



l/l«. 



< 4%. For example, while (001) is the softest direc- 
tion near Oeq and stays such upon biaxial expansion (Cu, 
Ni), it is one of the hardest in biaxially compressed met- 
als (Ag, Au, Cu) where (201) is the softest direction. 
Similarly, (111) is the hardest direction near the equi- 
librium and for <C Oeq, but it becomes softer than 

(110) and (201) in biaxially expanded Cu and Ni. Thus, 
there is a qualitative breakdown of the harmonic theory 
for strains of 4%, and presumably quantitative errors for 
even smaller strains. 

We also note similarities in the elastic bc^havior of these 
materials. Under expansion, both Cu and Ni exhibit 
strong softening of q{as, [001]) and somewhat weaker soft- 
ening of q{as, [111]), while q{as, [110]) becomes the elas- 
tically hardest direction. This order is reversed under 
biaxial compression of Ag, Au and Cu: q's for (001) and 

(111) harden, but the (110) and (201) directions soften. 



C. Discussion of anomalous softening of q{as,G) in 
terms of fcc/bcc energy differences 

The anomalous softening of q{as, [001]) in Ni and Cu 
for Us > Oeq reflects a small fcc/bcc energy difference for 
these materials. This can be seen by considering three 
energy surfaces that deform fee into bee: 

(i) £(7, y): The most general surface is the total en- 
ergy as a function of the tetragonal shear 7 and volume 



V, shown as contour in Fig. 4(a) for Cu. The tetragonal 
shear along (001) is defined by: 



(21) 



where c/o = y/2'yi . E{'-f, V) has (at least) three extremal 
points, denoted in Fig. 4(a) as solid dots: one correspond- 
ing to the fee state, one to the bee state and one to the 
bet state. These states obey the extremal conditions of 
vanishing derivatives: 




(22) 



Figure 4(a) shows that for Cu fee and bet are locally 

stable minima with respect to 7 and V, while bee is a 
saddle point (maximum with respect to 7 and minimum 
with respect to y).64,65 

(ii) Bain path E{pi): A more speciflc function £^(7) = 
^)ly=const deflned by the tetragonal Bain path,^^ 
connecting fee and bee structures. The Bain path is ob- 
tained by changing the c/a ratio while keeping V = ca^ 
constant. When c/a = 1 the lattice type is bcc and when 
c/a = \/2 it is fee. The energy as a function of 7 must 
have extremal points at both 7 values corresponding to 
the cubic symmetry fee (7fcc = 1) and bcc (7bcc = 2~3) 
states, as well as at least another bet point 7bct with 
a zero derivative £^'(7) = O.^^^'^^ Usually,'*^"''^'''^"^^ for 
fee stable elements the bcc lattice is unstable [i.e., E{'y) 



7 



has a local maximum at 7bcc] and the bet state (a local 
minimum) occurs for 7bct < 7bcc- 

(iii) Epitaxial Bain path E[ccq{as)]'- This deformation 
path is obtained by scanning c while Ug is kept fixed, 
which corresponds to epitaxial growth on a (001) sub- 
strate with lattice parameter Og. c is determined from 
the total energy minimization at a fixed a^: 

Eq.(23) defines the epitaxial path ^^(7), shown as a con- 
tinuous line in Fig. 4(a). Since c/og = ^2^^ and 
V — call A, this path implicitly relates the out-of-plane 
dimension c to the substrate lattice constant a^, much 
like Eq. (16) does in the harmonic case. As noted in 
Ref. 65, the epitaxial path crosses all extremal points 
of -E(7, V) because Eq. (23) is satisfied where condi- 
tions Eq. (22) hold. Therefore, if we parametrize the 
epitaxial strain energy along this path as a function of 
Os, it has a global minimum corresponding to fee, a lo- 
cally stable minimum corresponding to bet and a max- 
imum at the bee state, see Fig. 4(b). We see that as 
as increases from the equilibrium fee value, Cu sequen- 
tially passes through the bee and bet states where the 
strain energy Aii^''P'(as, [001]) is equal to the fee/bee and 
fcc/bct structural energy differences. When these en- 
ergy diflferences are much smaller than the characteris- 
tic values of the bulk deformation energies Ai?'^""^(as) 
[see Fig. 4(b)], then (/(os, [001]) is anomalously soft [since 
q{as, [001]) - ^E'"'^'^^y AE'^^^as) for a, = (2t/bcc)^]. 

In summary, the softness of q{as, [001]) for as > acq is 
a reflection of the geometric properties of the (001) epi- 
taxial deformation path (connection between cubic sym- 
metry fee and bee structures), and a small fee/bee energy 
difference, ASf='=/^<=<= < A£;^""'(a,). It is important that 
the fee and bee points correspond to lattices with cubic 
symmetry, since it ensures that the energy surface has 
extremal points there. In zincblende GaP and InP,^^ 
epitaxial (001) path has only one point of cubic symme- 
try (c/a = V2, corresponding to undistorted fee), and 
therefore the energy surface E{'j, V) is not required to 
possess additional extremal points. As a consequence, 
AiJ'^P' (as, [001]) is a monotonously increasing function of 
Us, and q{as, [001]) does not soften with increasing a^. 

The described mechanism also accounts for the soft- 
ening of q{as, [111]) for > a^q in Cu and Ni under 
biaxial (111) expansion, since this deformation takes fee 
{c/a = \/6) into bee (c/a = •\/6/4), albeit at a much 
larger strain. However, we have not found any simple 
structure corresponding to the compressive (201) strain 
which could explain the softening of q{as < agq, [201]) in 
Ag, Au and Cu. The latter seems to be caused by rel- 
atively loose packing of atoms within the (201) planes, 
imposing small energy penalty on decreasing the inter- 
atomic distances. Indeed, the nearest-neighbor distance 
in (201) plane is Og, compared to as/V2 in (111) or (001) 
planes with high values of q{as, G) for ag < a^q- 



IV. STABILITY OF SUPERLATTICES AND 
ALLOYS 

A. Constituent strain of super lattices 

The bulk formation enthalpy of superlattices [Eq. (8)] 
is expressed as a sum of the interfacial energy J(n, G) 
and constituent strain energy AE^{x,G). As given by 
Eq. (6), the latter is a weighted average of the epitaxial 
strain energies of coherently strained constituents, min- 
imized with respect to the common in-plane lattice vec- 
tors ai and sl2- For the high symmetry directions (001) 
and (111), these vectors are related by symmetry oper- 
ations of the superlattice, so that ai and a2 are propor- 
tional to the ideal fee unit vectors a? and via Eq. (2). 
Then Ai?^g(x, G) can be calculated by minimizing the 
following expression with respect to the superlattice pa- 
rameter asL: 

A£;°1(x,G) = min[(l - x)I\Ef{asi.,G) 

+ xAi;^P'(asL,G)]. (24) 

For lower symmetry directions G, the in-plane unit vec- 
tors ai and a2 may relax differently, and the angle 
C0S7 = ai-a2/|ai||a2| is also free to vary. For instance, in 
(110) superlattices, the vectors ai and a2 are not related 
by symmetry, and therefore may scale differently, i.e., in 
ideal fee |a°|/|a2| = -\/2 but in the superlattice generally 
|ai|/|a2| ^ \/2. Equation (24) is much simpler than the 
general Eq. (6) requiring minimization with respect to 
three degrees of freedom: lengths )ai|, )a2|, and the an- 
gle a = (ai,a2). In the present work we adopt Eq. (24) 
even for low symmetry directions, using the calculated 
AEf\as,G) from Sec. IIIB. 

Ai^Qg: Figure 5 shows the equilibrium constituent 

strain energies /^E^{x,G) for the size-mismatched Cu- 
Ag, Ni-Au and Cu-Au systems. They are determined 
from Eq. (24), using only the epitaxial strain energies 
AE'^g of the constituents. There are obvious similar- 
ities in AE^{x,G) for the three noble metal systems. 
(201) superlattices have the lowest constituent strain en- 
ergy below X w 0.2, after that (001) becomes the soft- 
est direction. (Ill) is the hardest direction over a wide 
composition range, except close to a; = 1 where (110) is 
slightly harder. 

This behavior can be explained by the properties of 
the epitaxial softening function q{as,G), discussed in 
Sec. IIIB. For example, consider Cu-Au from Fig. 5. 
Upon biaxial compression of Au (corresponding to x < 
0.5), q{as, [111]) increases rapidly (sec Fig. 2), increasing 
the elastic strain energy and making this an elastically 
hard direction. In contrast, q{as, [201]) for Au decreases 
with biaxial compression, and at x < 0.2 there is small 
energetic penalty for deforming Cu and Au to a common 
in-plane lattice constant. Increase of q{as, [110]) for Cu 
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with as eventually causes this to be the hardest direction 
in Au-rich Cu-Au superlattices. 

ashix)- Figure 6 shows the equilibrium in-plane lat- 
tice constant aQi^(x,G) that minimizes the constituent 
strain. These are also the equilibrium lattice param- 
eters for infinite period superlattices. The lattice pa- 
rameters asL(2^, G) show large deviations from Vcgard's 
law, with the behavior of osl(2;, [001]) being particularly 
anomalous. The very unusual composition dependence 
of the superlattice parameter for (001) deserves a closer 
scrutiny: At x w 0.2 the superlattice parameter changes 
discontinuously to the lattice parameter of the larger con- 
stituent. The constituent strain energy abruptly changes 
slope and settles down to a strictly linear composition 
dependence. Furthermore, Ai?^g(x, [001]) is very small 
in comparison with AiJ^g for other directions. These 
anomalies are direct consequences of the soft q{as, [001]) 
for biaxially expanded Cu and Ni, which in turn is a 
consequence of the small fcc/bcc and fcc/bct energy dif- 
ferences for these metals (Sec. IIIB). Indeed, for a suf- 
ficiently Au-rich system ^E"^^ is smaller than AE'^'^' fa- 
voring a superlattice constant close to the equilibrium 
lattice parameter of Au. This large lattice parameter 
happens to fall on the flat region of the strain energy 
curve around the bcc and bet states of biaxially expanded 
Cu (see Fig. 4), where a local bet minimum exists in the 
function on the right-hand side of Eq. (24), shifting down- 
ward in energy with increasing x. At some critical value 
of the composition, the minimum around uau becomes 
deeper than the minimum close to acu, which causes a 
discontinuous jump in asL- Loosely speaking, Cu de- 
forms all the way into the bet structure and Au does not 
deform at all. That also explains the linear decrease of 
A£'cs(a;, [001]) after the discontinuity, since AEf^ = 
and ASpPj = const in Eq. (24). 

In conclusion, we summarize the prerequisites for low 
elastic strain energy of infinite superlattices: 

(i) One of the components should exhibit a particularly 
soft elastic direction under biaxial strain, e.g., (001) in Cu 
upon epitaxial expansion and (201) in Au upon biaxial 
compression. 

(ii) The lattice mismatch between the constituents 
should be large enough to explore the regions of anoma- 
lous softness. 

We stress that the unusual behavior shown in Figs. 6 
and 5 (crossing of different directions, discontinuities, dif- 
ferent skewnesses of AE^(x,G) curves) are caused by 
the anharmonic q{as,G), and cannot be obtained within 
the harmonic theory with lattice parameter independent 

9harm(G).^^ 



B. Describing chemical interactions via the 
mixed-space cluster expansion 

The energy of a bulk alloy AH^^^^{x) of Eq. (5), and of 
an epitaxial alloy Aff^^(x) of Eq. (11) cannot be com- 
puted directly from LDA since configurationally random 
structures are involved. The approximate approach is 
either large supercells or a first-principles mixed-space 
cluster expansion. In the latter approach, a spin vari- 
able Si is assigned to each lattice site Ri which takes a 
value +1 if the site is occupied by an atom of type A, or 
— 1 if the site is occupied by an atom of type B. The for- 
mation enthalpy of an arbitrary structure a is expressed 
in the following form: 



AJfcE(a) =^Jpair(k) \S{k,a)f 



MB 

+ J2DfJfIifi<j)+AEcs{a). (25) 
/ 

where J(k) is the Fourier transform of real-space pair 
interactions and S'(k, a) is the structure factor, 

Jpair (k) = ^ Jpair (Ri - R,' )e-^''^^ , (26) 
j 

S{k,a) = J2Sje-''^'''. (27) 
j 

The second sum in Eq. (25) nms over symmetry in- 
equivalent clusters constituted by three or more lattice 
sites. Df is the number of equivalent clusters per lat- 
tice site, and 11/(0-) are structure-dependent geometrical 
coefficients (simple lattice averages of the cluster spin 
products). The last term in Eq. (25) is the constituent 
strain energy AEcs{o') oi the structure a. It is designed 
to reproduce the elastic strain energy of coherent long- 
period superlattices^® which are calculated directly from 
the LDA (see Sec. IV A): 

AEcs{(T) = J2Jcs{x,k)\S{k,a)\\ (28) 

k 

, , AEpJx,k) 

The equilibrium constituent strain energies AE^{x,k) 

have been deduced from the directly calculated 
AE'^P^{asL,G) minimizing Eq. (24) with respect to the 
common in-plane lattice constant asL- They are fitted 
by series of Kubic harmonics with composition depen- 
dent coefficients ci{x): 



AEcs{x,G) = Y,ci{x)Ki{G), 



(30) 



which are used to evaluate AEq^ [x, G) for any direction 
G. Structure factors S'(k, a) in the long-period superlat- 
tice limit are nonzero only for k — > 0, where Jcs{x,k) 
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is a nonanalytic function of k, reflecting the directional 
dependence of the constituent strain energy. 

The effective cluster interactions Jj and Jpair(k) are 
determined by fitting Eq. (25) to a large number (30 
to 40) fully relaxed first-principles LDA formation en- 
thalpies of simple ordered structures. Most of these 
ordered structures are short-period superlattices along 
(001), (111), (110), (201) and (113).*^° The calculations of 
T = total energies employ the full-potential linearized 
augmented plane wave method^^ (FLAPW) and local 
density approximation (LDA) for the electronic exchange 
and correlation. The total energy is minimized with re- 
spect to all structural degress of freedom, i.e. both the 
atomic positions and cell-external coordinates are fully 
relaxed. Complete discussion of the LDA calculations 
and cluster expansions for Ag-Au, Cu-Ag, Cu-Au and 
Ni-Au can found in Ref. 60. 



C. Stability of finite period metal superlattices 

Having obtained all ingredients of Ai?cE(o') [Eq. (25)] 
from LDA calculations on small unit cell structures, we 
can use this equation to predict the energy of any con- 
figuration <T, in particular superlattices. Figure 7 shows 
the bulk formation energies of (j4)„/(i?)„ superlattices 
for the studied noble metal systems. The interfacial en- 
ergies I{n,G), extracted from AHsL{n,G) according to 
Eq. (8), are shown in Fig. 8. Several interesting observa- 
tions can be made from these figures: 

(i) I{n, G) are approximately constant after n > 5. 

(ii) For ordering systems (Cu-Au and Ag-Au), the in- 
terfacial energies are negative (see Fig. 8). Negative in- 
terfacial energies are the cause for the upward slope of 
AHsh{n,G) curves in Fig. 7. This indicates a chemical 
preference for having unlike atoms at the interface and 
a tendency to form ordered structures at T = 0. For 
instance, LIq, the observed ground state of CuAu, is a 
monolayer (Cu)/(Au) superlattice along (001). The or- 
der of A7?sL(n,G) is lowest (001) and highest (111) for 
Cu-Au, and lowest (110) and highest (111) for Ag-Au 
superlattices. 

(iii) For the phase separating Cu-Ag, all interfacial en- 
ergies are positive. AHsh{n,G) decreases with n for all 
directions and reflect the tendency to coherent phase sep- 
aration over ordered superlattice formation. Interfaces 
between Cu and Ag are energetically very costly. The 
order of AHsL{n,G) is again lowest (001) and highest 
(111). 

(iv) Ni-Au has the most interesting behavior of 
AHshirijG) and I{n,G). It exibits phase-separating 
type AHsh{n, [001]) (decreasing with n), ordering type 
AHsL{n, [110]) (increasing with n), and a nearly con- 
stant AHsi,{n, [111]). Does this mean that interfaces 
in some directions are energetically favorable, while in 
other directions they are energetically costly? The an- 
swer is: No. In Ni-Au, just like in Cu-Ag, all isolated 
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interfaces have positive formation energies. However, the 
interaction between the interfaces along (110) is strongly 
attractive in Ni-Au, and leads to a net chemical energy 
gain for some short-period superlattices. Indeed, Fig. 8 
shows that all interfacial energies of Ni-Au are positive 
in the limit n — !■ oo (when there is no interaction be- 
tween the interfaces), but decrease for short periods and 
are negative for (110) n < 3 superlattices. As we show 
in Ref. 61 the competition between the constituent strain 
energy, interfacial energy I{n ^ oo, G) and ordering- type 
interaction between the interfaces leads to unusual short- 
range order in Ni-Au solid solutions. 

(v) It is interesting that in the phase separating Ni-Au 
and Cu-Ag the lowest interfacial energy 7(n oo,G) oc- 
curs for the close-packed {111} interfaces, and the highest 
for {001} interfaces. This situation is completely differ- 
ent in the ordering systems Cu-Au and Ag-Au, which 
have {110} as the lowest and either {111} or {001} as 
the highest I{n oo, G). 

(vi) Figure 7 shows the enrgies of the random alloys at 
the equiatomic composition. We see that in Cu-Au and 
Ag-Au all long-period superlattices are unstable with re- 
spect to the formation of a random alloy. In Ni-Au the 
random alloy is less favorable than coherent phase sepa- 
ration in the (001) direction, but slightly more favorable 
than infinite coherent superlattices along (110) and (111). 
However, short-period (110) superlattices are lower in en- 
ergy than the random alloy. All (111) superlattices of 
NiAu have higher formation enthalpies than the random 
alloy. In Cu-Ag only the long-period (001) superlattices 
have lower bulk formation enthalpies than the random 
alloy. The epitaxial growth of CuAg and NiAu alloys is 
discussed more thoroughly in Sec. IV E. 

(vii) In size-mismatched systems (Cu-Ag, Cu-Au, and 
Ni-Au) AHshinjC) exhibit the same order as the con- 
stituent strain AE^{x, G), i.e., AHsj,{n, [001]) is lowest 
and AHsL,{n, [111]) is highest. It suggests that low con- 
stituent strain stabilizes even short-period superlattices. 
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FIG. 8. Interfacial energies of Cu-Au, Cu-Ag, Ni-Au and 
Ag-Au. 



D. Comparison of the trends in stability of metal 
and semiconductor superlattices 

Growth of semiconductor superlattices is a more ma- 
ture are than than growth of metal superlattices, and 
much more data are available at present. Thus, it is of 
interest our results in Figs. 7 and 8 for metals with anal- 
ogous results for semiconductors. ^^'^^ 

Lattice-mismatched semiconductors generally have 
AH^^{x) > and AH^l^^ > 0. Thus, they re- 
semble Ni-Au and Cu-Ag rather than the compound- 
forming system Cu-Au. LDA calculations reveal that 
AH^l^^{n,d) for G ^ (111) and G = (001) look exactly 
like in Cu-Ag or Ni-Au: the energy decreases as the pe- 
riod n increases, and the interfacial energies are mostly 
positive. However, in the (110) and (201) directions, 
the interfacial energies are negative, and Ai7gL^'^(n, G) 
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FIG. 9. Mixing enthalpies A/fmix (in meV/atom) for bulk and epitaxial equiatomic Cu-Ag and Ni-Au alloys. All epitaxial 
calculations assume that the substrate is lattice matched to the random alloy. AEcs is the sum of epitaxial strain energies of 
pure elements [see Eq. (11)]. 



increases with n, like in Ni-Au and Cu-Au. Hence, 

semiconductor supcrlattices behave generically as Ni-Au 
superlattices. However, short-period (201) semiconduc- 
tor superlattices (e.g., the chalcopyrite structure, corre- 
sponding to n = 2) have a lower energy than the random 
alloy, while in Ni-Au it is the (001) short-period super- 
lattices that have lower energies than the random alloy. 
Hence, while the Ni-Au random alloy can lower its energy 
by developing (001) ordering, semiconductor random al- 
loys can lower their energy by developing (201) ordering. 
Both in Ni-Au and semiconductor alloys, the ultimate 
ground state is incoherent phase separation. 

E. Epitaxial growth and surface intermixing 

Recent experimental studies"^^'"^^ have grown epitax- 
ial films of noble metal alloys which are immiscible in 
the bulk form. For instance, Stevens and Hwang''^ have 
grown Cu-Ag alloys on a Ru(OOOl) substrate, demon- 
strating that Cu and Ag are miscible at T = 823 K, 
despite the fact that in bulk, Cu and Ag are strongly 
immiscible at this temperature and composition. It has 
also been observed that Au deposited on Ni(llO) sur- 
face replaces it in the first surface layer forming a sur- 
face Ni-Au alloy,^^ although Au is completely insoluble 
in bulk Ni. In what follows we show that the stabiliza- 
tion of epitaxial solid solutions from bulk- immiscible con- 



stituents can be explained by the additional destaUKza- 
tion of the constituents due to the epitaxial constraint. 
Indeed, Eq. (11) shows that the epitaxial mixing enthalpy 
"^-^mix considerably lower than the bulk mixing 

enthalpy AH^^^ if the sum of the constituent strain en- 
ergies on the right hand side is large. 

Figure 9 shows the results for the epitaxial stabilization 
of equiatomic NiAu and CuAg alloys, assuming that the 
substrate is lattice matched to the disordered alloy. 

(i) Disordered CuAg and NiAu alloys have large pos- 
itive bulk mixing enthalpies AH^^, in agreement with 
the observed bulk imniiscibility. 

(ii) Epitaxy destabilizes the constituents, and hence 
stabilizes the epitaxial alloy in all cases. This effect is 
much larger for the elastically hard direction (111) than 
for the soft (001) direction. 

(iii) The epitaxial mixing enthalpy SH^^ for (111) be- 
comes negative in Ni-Au, showing that the solid solution 
is energetically favored over the epitaxially phase sepa- 
rated state. In CuAg, ^^111 positive and these 
alloys are unstable under epitaxial conditions at T = K. 

(iv) Epitaxial conditions lead to a significantly en- 
hanced miscibility since SH^^ < SH^^. A simple 
mean-field estimate of the miscibility gap temperature for 
CuAg grown on a nearly lattice-matched Ru(OOOl) sub- 
strate [equivalent to a fcc(lll) substrate] gives Tmg = 
2AH^}^ = 150 K. Thus, for (lll)-epitaxy at the tem- 
perature (823 K) of Steven's and Hwang's experiment. 
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our calculations predict complete solubility of Cu-Ag, as 

observed. 

(v) The epitaxial stabilization is strongly dependent on 
the substrate orientation. A bigger effect can be observed 
for elastically hard directions, e.g., (Ill) and (110) for 
noble metal alloys. 

V. SUMMARY 

We have investigated the effects of anharmonic strain 
on the stability of epitaxial films, superlattices and epi- 
taxially grown disordered alloys. We find that anhar- 
monic epitaxial strain produces certain qualitative and 
quantitative features absent in the harmonic theory. In 
particular, 

(i) Epitaxial softening functions q{as , G) are strongly 
dependent on the substrate lattice constant ag, while 
they are constants in the harmonic theory. For instance, 
as a consequence of the small fcc/bcc and fcc/bct en- 
ergy difference, biaxially expanded Cu and Ni show dras- 
tic softening of q{as, [001]). Furthermore, biaxially com- 
pressed Cu, Ag, and Au have low values of g(as, G) along 
directions (201) and (110) with relatively loose packing 
of atoms in the epitaxial planes. 

(ii) The dependence of (/(a^, G) on the direction G can 
differ from harmonic predictions. For instance, (110) is 
the hardest direction in biaxially expanded Cu and Ni, 
and (201) is the softest in biaxially compressed Cu, Ag 
and Au. The harmonic formula always predicts either 
(111) as the hardest and (001) as the softest direction, 
or vice versa. 

(iii) The strain energy of infinite coherent superlat- 
tices exhibits marked anomalies associated with the an- 
harmonic behavior of constituent q{as,G). The size- 
mismatched systems Cu-Ag, Cu-Au and Ni-Au exhibit 
very low constituent strain for Ag- and Au-rich (001) su- 
perlattices, since (001) is the easy direction for biaxial 
expansion of Cu and Ni. Similarly, (201) superlattices 
with small Ag or Au content have low coherency strain 
energies because this is the easy deformation direction for 
biaxially compressed Ag and Au. The in-plane lattice pa- 
rameter asL of long-period (001) superlattices suffers a 
discontinuous jump around x ~ 0.2, and other directions 
show considerable deviations from linear behaviour. 

(iv) These elastic anomalies are less pronounced in 
short-period superlattices, although they contribute to 
the structural stability of (001) superlattices. Short- 
period bulk superlattices are stable in Ag-Au and Cu-Au 
due to negative interfacial energies. Ag-Au and Ni-Au 
have positive interfacial energies, leading to superlattice 
formation being energetically unfavorable with respect 
to phase separation. The interaction energy between in- 
terfaces in Ni-Au is so strong that short-period (n cx 2) 
superlattices along (110) are more stable than the long- 
period superlattices with fewer interfaces. 



(v) Epitaxially grown disordered alloys can be stabi- 
lized even if the system phase separates in bulk form. 
This effect is caused by additional destabilization of the 
phase separated state due to the epitaxial constraint on 
the constituents, requiring them to be coherent with the 
substrate. The stabilization is more pronounced for elas- 
tically hard directions with high values of q{as,G), e.g. 
(111). For instance, we find that even though Ni-Au 
and Cu-Ag phase separate in the bulk {AH^}^{x) > 0), 
equiatomic Nio.sAug.s alloys are miscible when grown on 
a lattice-matched (111) substrate, while Cuo.sAgo.s on a 
(111) substrate is immiscible at T = K but miscible at 
T > 150 K. Neither Nio.sAuo.s nor Cuo.sAgo.s are mis- 
cible when grown on a lattice-matched (001) substrate, 
corresponding to a low energy penalty on the phase sep- 
arated constituents. 
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